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ABSTRACT 

Aims. Coronal seismology is extensively used to estimate properties of the corona, e.g. the coronal magnetic field strength are derived 
from oscillations observed in coronal loops. We present a three-dimensional coronal simulation including a realistic energy balance in 
which we observe oscillations of a loop in synthesised coronal emission. We use these results to test the inversions based on coronal 
seismology. 

Methods. From the simulation of the corona above an active region we synthesise extreme ultraviolet (EUV) emission from the 
model corona. From this we derive maps of line intensity and Doppler shift providing synthetic data in the same format as obtained 
from observations. We fit the (Doppler) oscillation of the loop in the same fashion as done for observations to derive the oscillation 
period and damping time. 

Results. The loop oscillation seen in our model is similar to imaging and spectroscopic observations of the Sun. The velocity 
disturbance of the kink oscillation shows an oscillation period of 52.5 s and a damping time of 125 s, both being consistent with 
the ranges of periods and damping times found in observation. Using standard coronal seismology techniques, we find an average 
magnetic field strength of B kink =79 G for our loop in the simulation, while in the loop the field strength drops from some 300 G at 
the coronal base to 50 G at the apex. Using the data from our simulation we can infer what the average magnetic field derived from 
coronal seismology actually means. It is close to the magnetic field strength in a constant cross-section flux tube that would give the 
same wave travel time through the loop. 

Conclusions. Our model produced not only a realistic looking loop-dominated corona, but also provides realistic information on the 
oscillation properties that can be used to calibrate and better understand the result from coronal seismology. 

Key words. Sun: corona — Sun: activity — Sun: UV radiation — Sun: oscillations — Magnetic fields — Magnetohydrodynamics 
(MHD) 


1. Introduction 


Loops dominate the appearance of the corona of the Sun. 
In particular in active regions observed at extreme ultravio¬ 
let (EUV) wavelengths they are seen as fine threads outlining 
the magnetic field. In response to a significant localized en¬ 
ergy de position, like in a flare, these loops can bee seen to os¬ 
cillate (INakariakovet al.l l l999 ). Since then oscillations in the 
corona have been studied extensively in terms of theoretical in¬ 
ves tigations, numerical models, and observations (e.g. reviews 
by Nakar iakov & Verwicht ell2005[ [Baneriee et alJl2007t IWang 
ariakovll^ 


201 It iDe Moortel & Nakariakovl 1201 21) . These oscillations of 
coronal loops are used for diagnostics of corona loops, in par- 
tic ular of the magn etic field strength. This w as proposed already 
by lUchidal (119701) and iRoberts et alJ (119841) who provided the 
theoretical ground work. Traditionally EUV imaging and spec¬ 
troscopy gave access only to plasma properties, i.e., temperature, 
density, abundance and flows (e.g.[Maris kalfT992h . Coronal seis¬ 
mology holds the key to infer information also on the magnetic 
properties, in particular for the field strength and resistivity. The 
key techniques for coronal seismology give access to the mode 
of the oscillation and the phase speed of the corresponding wave. 

The oscillations of coronal loops have been investigated 
mainly by two techniques. The first techniqu e measu res the dis- 
place me nts and disturbanc es in EU V images ( AschwandgnetajJ 


19991: I Nakaria kov et al.l 
Nakariakov & Ofmanl 12001 


19991 : IDe Moortel et al. 2000: 


Aschwanden & Schriiven 


2011 


Guo et al 


Yuan & Nakariakovl 

20151) . 


1201 2t IVerwichte et al] 1201 3bt 


riodic 


The second technique investigates the pe- 
in the Doppler shifts obtained from EUV 


spectrometers (Ofman & Wang) 2002; IWan^ et al. 

20031 

2007; Van Doorsselaere et al. 2008; Erdelyi & Tarovan 

2008; 

Ofman & Wang 2008; 

Mariska & Muglach 201C4 Tian et al. 

2012|). Both techniques give an average magnetic field strength 


of typically 10 G to 100 G in the loop. This is in general 
consistent with the field strength derived fro m extrap olations 
of the photospheric magnetic field (e.g. ISchriiver et alJ 120061 : 
IWiegelmann & Sakuraii 12012 ^. However, the results from 
coronal seismology can only provide some average value of 
the magnetic field along the loop, while clearly the magnetic 
field has to expand with height and therefore the magnetic field 
strength will change along the loop. Consequently, it remains 
unclear what this average derived from coronal seismology 
really means. 


While the original work was assuming a rather simple setup 
(e.g. IRoberts et al.ll984l) . more recent theoretical efforts have ac¬ 
counted also for the more complex structure of the real Sun, e.g., 
curved geometry, density stratification , or non-uniform cross 


section (|Van Doorsse 

aere et al. 2004; Verwichte et al. 2006; 

Erdelyi & Verth 20071 

Arregui et al. 2007'; Goossens et al. 2009; 

Ruderman & Erdelyi 2009: Selwa et al. 2011). In a 3D model 

De Moortel & Pascoe 

(2009!) investigated the estimate of the 


magnetic field strength by coronal seismology. In their model 
the magnetic field strength and number density along the coro- 
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nal loop are assumed to be constant dPascoe et al.ll2009 ). This 
configuration sets a single reference value, and the authors did 
find a difference between their reference value and field strength 
derived from the oscillation of the model structure. 

To make things worse, in the real corona the magnetic field 
strength varies along the loop, as does the number density, de¬ 
spite of the large density scale-height in the corona, and this 
migh t have a significant impact on the results of coro n al seis mol¬ 
ogy dOfman et al.ll2012l) . lAschwanden & Schriiverl (1201 lh and 
IVerwichte et al.l (12013ah compared the magnetic field strength 
from coronal seismology with those obtained from potential or 
force-free extrapolations of the magnetic field in the same struc¬ 
tures. Because of the limitation of the approaches (i.e. poten¬ 
tial or force-free), it is uncertain if the extrapolated magnetic 
field line actually match the obser ved loops. Promising exam¬ 
ples were presented by iFeng et ail ( 2007 ) who showed that the 
magnetic field lines in a linear-force-free extrapolation do follow 
the loop structures reconstructed from stereoscopic observations. 

Given all the difficulties of getting the reference values from 
observations, it is highly desirable to test coronal seismology 
in a model corona, which has realistic plasma properties and a 
magnetic field configuration similar to tha)_of a real active re- 
gion. Forward corona l models dGudiksen & Nordlundll2005allbl : 
IBingert & Peter! 1201 lh account for the cooling through the op¬ 
tically thin radiation and the highly anisotropic heat conduction 
along magnetic field lines, and su cces sfully reproduce a vari¬ 
ety of coronal phenomena dPeterl 120151) . In those models the 
the horizontal motions in the photosphere ind uce currents in 
the corona, either through fieldline braiding (iParkerl R 972l) or 
flux-tube tectonics (iPriest et al.ll2002l) . The Ohmic dissipation of 
these currents is sufficient to heat the coronal plasma to over 
one million K. The energy distribution in this type of model 
is co n siste nt with the expectation of the nanoflare mechanism 
(IBingert & Peterl201 3i) . The full treatment of the energy balance 
allows these models to resemble the plasma properties in the real 
corona, so that the synthesised emission from these models can 
be directly compared with real observations. These model suc¬ 
cessfully explain some basic f eatures o f coronal loops (e.g. the 
non-expanding cross section, iPeter & Bingertll2012 1 A one-to- 
one data-driven simulation can reproduce the appearance and 
dynamics in the partic ular solar active region that drives the sim¬ 
ulation (iBourdin et al.ll2013h . 

Being successful in modelling the plasma properties and 
general dynamics of the corona, a further challenge is how 
well the loop oscillation in a realistic model resembles that in 
real observations. We analyse the coronal emission synthesised 
from a 3D MHD model in which we clearly see loop oscil¬ 
lations. Analysing the synthetic data in the same way as ob¬ 
served coronal oscillations, we estimate the average magnetic 
field strength in the loop. Having access to the full 3D cube 
of the simulation data, we can compare this average value to 
the actual magnetic field strength that varies along the model 
loop. This comparis on is a complementary to previous nu- 
meric al experiments ( DeMoorteL&Pascoe 2009) and observa¬ 


tions (lAschwanden & Schriiveil 120111 : IVerwichte et alJ l2013al) . 
and helps to better understand the implications derived from the 
field strength inferred from coronal seismology. 


2. Model setup 

The coronal model we analyse here is based o n the mod¬ 
elling strategy as described by IBingert & Peteii (1201 ll) . The 
m odel setup i s the same as in our previous model described 
in IChen et al. ( 2014 1 1201^ ). However, the simulation described 


here has significantly higher resolution: The 147 x 74 x 50Mm 3 
volume is now resolved by 1024 x 512 x 256 grid points, which is 
an increase by a factor of 4 in each horizontal direction. The grid 
spacing is uniform in the horizontal direction (144 km grid spac¬ 
ing) and non-uniform in the vertical direction (smoothly chang¬ 
ing from 30 km in the photosphere to 190 km in the coronal part, 
ranging from 2 Mm to 4 0 Mm). To so lve the full MHD equations 
we use the Pencil Code (iBrandenburg & Dobleii2002h n 

Our model corona is driven by the emerging magnetic field 
and the flows in the photosphere. These are taken from a simu¬ 
lation of the emergence of a magnetic flux tube from th e uppe r 
convection zone through the photosphere (ICheung et alJl201H) . 
Here we use a model where the emer ging flux tube has no 
imposed twist (iRemnel & Cheun5l2014l) . In the process of the 
emergence a pair of sunspots is formed by coalescence of small 
magnetic patches. This flux emergence model covers only a 
small part of the the photosphere. We take the time-dependent 
output of the flux emergence simulation at the solar surface 
(magnetic field, velocity, density and temperature) and impose 
thi s at the lower boundary of our coronal model (just as in 
IChen et aPl20l4 but now at higher resolution). 

The major limitation for the time step in the explicit time¬ 
stepping scheme of our model is due to the Spitzer heat conduc¬ 
tion. Therefore we evolve the equations in an operator-splitting 
manner. The governing equations without the heat conduction 
term are evolved by a regular 3rd-order Runge-Kutta scheme. 
We then evolve an equation for the the heat conduction term 
alone in a sub-cycle. In the sub-cy cle, we use a super-time- 
stepping scheme dMever et al.ll2012l) . which allows us to evolve 
the energy equation with a time step much larger than that for 
the original Runge-Kutta scheme. 

The second-smallest time step in the simulation is (mostly) 
the Alfven time step. The Alfven speed in the model corona can 
be very large (up to 15 000km/s), which would then limit the 
time step. Fortunately these large speeds are found only in a 
small fraction of the computational domain. To overcome this, 
we control the Alfven speed v A b y limit ing the Lorentz force 
in the same way as IRemnel et al.l d2009l) . By multiplying the 
Lorentz force in the momentum equation by a correction fac¬ 
tor /a we ensure that the resulting effective Alfven speed in the 
model, 

VA = fl /2 V K , ( 1 ) 

is al ways smaller than a maximum speed u rna x=2000kms _1 . As 
in lRempel et al.l d2009l ) the correction factor is defined as 

/a = ^maxK + ^max) ■ ( 2 ) 

For a typical coronal loop with a number density of 10 9 cm -3 
and a magnetic field strength of 50 G at its apex, the Alfven 
speed («3000kms -1 ) is reduced by about a factor of 1.5. Sig¬ 
nificant corrections will mostly occur in low-density regions 
(where the coronal Alfven speed gets large), which will not ap¬ 
pear bright in the synthetic images and thus can be expected to 
play a minor role only for the analysis we present in this paper. 
Despite this limiting of the Alfven speed, in the model corona 
the plasma-/? will still be well below 1 and the effective Alfven 
speed will still be an order of magnitde larger than the maximum 
sound speed (250km s _1 at 3MK). In summary, this limiting of 
the Alfven speed will have a minor effect only on our results, 
but it will provide a significant speed-up of the numerical simu¬ 
lation. 

1 See also http://pencil-code.nordita.org. 
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Together, the treatment of the heat conduction and the Alfven 
speed, introduced only to speed-up the simulation time, provide 
a speed-up of typically a factor of five, at little extra computa¬ 
tional cost. 

Our high-resolution simulation is able to resolve small-scale 
photospheric magnetic structures and flows related to the granu¬ 
lar motions in the photosphere which are used as input from the 
flux-emergence simulation. The interaction of the plasma flow 
and magnetic field in the model photosphere, especially at the 
outer edge of strong flux concentrations, produces the enhanced 
upward Poynting flux that can bring enough energy into the up¬ 
per solar atmosphere and power the coronal loops. This energy 
flux sustains a more than 3 MK hot corona at pressures according 
to the classical scaling laws dRosner et al .11 19781) . 

We show in Fig.Q] a snapshot of a time series of the syn¬ 
thesised EUV images according to the response function of the 
211 A channel of the Atmospheric Imaging Assembly (AIA, 
iBoerner et al.l 20121) . Here we follow the procedure outlined in 
dPeter & Bingerdl2012l) to calculate the AIA emission expected 
from the model corona. In these synthetic images we find nu¬ 
merous bright EUV loops, which represent the coronal plasma 
at temperatures around 2 MK. In these loops the number density 
is about K) 1 'cm 3 . The properties of the coronal plasma are in 
quantitative agreement with our previous simulation, and they 
are also consistent with typical values derived from real obser¬ 
vations. 

Beyond the appearance in a synthetic EUV snapshot, the 
model also captures the dynamic nature of the real corona to 
some extent. The animation associated with Fig.Q] shows that 
bright features can show up or disappear within a few minutes in 
a certain EUV passband. This is consistent with modern obser¬ 
vations, for example from AIA. 


3. Oscillation observed in synthesised coronal 
emission 

3.1. Trigger of the oscillation 

Similar as in our previo us model of an emerging active region 
(IChen et all 1201 4 1201 5l) loops form in the active region in re¬ 
sponse to the heating driven by footpoint motions in the pho¬ 
tosphere. This leads to a more or less continuous evolution of 
the loops: some loops form, others fade away. A snapshot of the 
loops in coronal emission is displayed in Fig.Q] the temporal 
evolution is shown in the attached movie. In the following we 
relate all times relative to the time 72.5 minutes after the actual 
start of the simulation. 

At the time fss200 s some loops start to oscillate just after a 
strong transient brightening of some close-by loops. These lo¬ 
cations are indicated in the snapshot in Fig.Q] as well as in the 
attached movie. 

We traced the cause of the strong transient brightening as 
being due to an enhancement of Ohmic heating along the field 
lines at the flank of the active region. The cool plasma is quickly 
heated to a temperature at which the response function of the 
AIA 211 A channel peaks. The increase of the heating is caused 
by an increase of the currents in response to reconfiguration of 
the coronal magnetic field driven by photospheric flows in the 
periphery of the sunspots (as the result of the flux emergence). 
The strong brightening associated with this increased heating 
rate can be considered as a very small flare on the Sun (even 
though we do not imply here that this is a flare model). 

The increased currents come along with a transient increase 
of the Lorentz force, leading to a kick perpendicular to the 



Fig. 1 . System of coronal loops synthesised from 3D MHD model. 

This displays a snapshot of the model corona as it would appear in 
an EUV image taken by AIA in the 211 A channel (in logarithmic 
scale). It is dominated by emission from Fe XIII showing plasma at 
around 2MK. The distance between the two footpoints of the loop 
system is about 35 Mm and the loops have lengths of about 45 Mm to 
50 Mm. The arrows indicate the position of the brightening that triggers 
the oscillation (see Sect. l3.lb and the oscillating loop (Sect. l3.3V The 
full temporal evolution over 54 minutes is available as a movie in the 
online edition. The movie starts early in the simulation, when there is 
no coronal emission. In response to the flux emergence coronal loops 
form, and at about 72.5 minutes a trigger sets the oscillation in motion. 

Then a second counter in the movie shows the time in s starting at 72.5 
minutes, consistent with the time used in Fig.Q] 

The movie is also available at 

http://www2.mps.mpg.de/data/outgoing/peter/papers/2015-osci/movie-figl.mp4 


field lines. The disturbance propagates from the brightening site 
across the active region and triggers the transverse oscillation 
of the nearby loops. Again, this can be considered in analogy 
to eruptions on the real Sun, which are suggested to be the main 
trigger mechanism of observed loop oscillations. This oscillation 
is clear, albeit not violent, in the movie attached to Fig.Q] 

The travel time of the triggering disturbance from the tran¬ 
sient brightening to the farthest visible loop is about 12 s (prop¬ 
agation with an Alfven speed of almost 2000kms _1 across a 
distance of some 25 Mm). This is why at the cadence of 10 s for 
the movie shown with Fig.Q]the oscillation of the loop seems to 
start at almost the same time as the transient brightening trigger¬ 
ing the event. 

The actual process triggering the oscillation is a very inter¬ 
esting process in itself, but it is not the main interest of this study. 
In contrast, we concentrate here on the observable consequences 
of this event and investigate to what extent the methods of coro¬ 
nal seismology can recover the atmospheric conditions at the lo¬ 
cation of the oscillating loop, in particular the magnetic field. 

3.2. imaging and spectroscopy of model data 

The synthesised imaging data (movie with Fig.Q) reveal a trans¬ 
verse oscillation the loop that lasts for about 300 s and de¬ 
cays gradually. Its appearance is similar to the oscillations 
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widely found in EUV observations. Mostly this type of os¬ 
cillation _is interpreted as the standing, fast kink mode (e.g. 
iNakariakov & Verwichtel 120051) . Although the loop oscillation 
is clear in the animation for the human observer, a quantitative 
analysis of this oscillatory displacement of the loop position is 
not so clear, because the amplitude of the oscillation is only of 
the order of or even smaller than the width of the coronal loop. 

Alternatively, one can investigate the velocity disturbances 
through Doppler shifts from (synthesised) spectroscopic obser¬ 
vations. Because the loops are inclined, when observed from 
straight above the transversely oscillating loops will result in 
periodic Doppler shifts. To study this, we simulate an observa¬ 
tion from the top of the active region, i.e. as if we would ob¬ 
serve an active region at the disk center. For this we calculate 
the emission line profile at each grid point based on the output 
of the MHD simulation and then inte grate along a vertical line- 
of-sight. This procedure follows [Peter et al.l d2004l2006l) . To be 
consistent with the imaging data shown in Fig. [I] we study line 
profiles from the same ion and choose the Fe XIII line at 202 A 
line that has been o bser ved abu ndantly with the EUV Imaging 
Spectrometer tEISi lCulhane et akll2007h . 

The Fe XIII profiles are mostly close to a single Gaussian. 
Their width is determined by the local plasma temperature and 
the distribution of the line-of-sight velocities. To derive the line 
intensity and shift we take the zeroth and first moment of the line 
profile. The resulting Doppler map in Fig.[2]basically shows the 
vertical (line-of-sight) velocity of the 2 MK hot plasma in the 
active region (at the same time as the snapshot in Fig.QJ. This 
Doppler map shows the line shifts at the same time across the 
whole map. In a real observation with a slit spectrometer one 
would have to produce a raster scan to obtain this map, e.g. EIS 
would need up to one hour to raster a field-of-view as shown in 
Fig.[2](depending on observation parameters). Thus a real obser¬ 
vation would look quite different than this instantaneous Doppler 
map. 

In a real observation usually one would perform a so-called 
sit-and-stare observation to catch the oscillation, i.e., one would 
keep the slit at a fixed position in the active region and study the 
temporal evolution at that location. For our synthetic observa¬ 
tions, we choose a slit oriented along the ^-direction located just 
in the middle between the footpoints of the loop (solid line in 
Fig.0. The dotted line indicates the location of the loop that is 
seen to be oscillating in the intensity data (cf. movie with Fig.|T}. 
Here we catch the loop in a phase of the oscillation moving away 
from the virtual observer. The slit roughly crosses the apex of the 
loops. 

For the further analysis we extract the line intensity and the 
Doppler shift along this slit from the synthesised spectral data 
with a cadence of 10 s. This roughly matches the typical cadence 
used in real observations. The resulting time series of this syn¬ 
thetic observation is comparable to real data acquired from active 
region loops in a sit-and-stare observation. 


3.3. Quantitative analysis of the loop oscillation 

To analyse the oscillation we first check time-space diagrams 
for the line intensity and Doppler shift. Just like fo r sit- 
and-stare observations of the real Sun (e.g. IWang et ail 120091 : 
iMariska & Muglachl2010l) these show the line intensity and shift 
as a function of time and space (along the slit). For the slit posi¬ 
tion indicated in Fig. [2] these are displayed in panels (a) and (b) 

of Fig. [3] 



40 60 80 100 

x [Mm] 


Fig. 2. Map of Doppler shifts of synthesised Fe XIII 202 A data. This 
shows the active region seen from straight above, i.e. the line-of-sight 
is vertical, corresponding to an observation at disk center. The black 
solid line indicates the position of the slit used to simulate a sit-and- 
stare observation for the analysis of the Doppler oscillation (see Fig.[3] 
and Sect. l3.2b . The dotted line indicates the location of the oscillating 
loop as seen in the movie attached to Fig.Q] 


The line intensity does not show a clear oscillation in the y- 
direction in this diagram. This is consistent with the impression 
from the synthetic AIA images that the amplitude of the oscil¬ 
lation (in space) is smaller than the loop width. In addition, the 
loop is inclined and the oscillation is transverse, which further 
reduces the amplitude of loop displacement. 

In contrast, the oscillation is very clear in Doppler shift 
(Fig.Eb). It starts at around £=200 s, which is consistent with 
the oscillation seen as a slight displacement of the loop in coro¬ 
nal emission (see movie with Fig.|T}. The oscillatory pattern is 
seen over the whole field-of-view shown in Fig.JTjb underlining 
the impression from the intensity images that the trigger leads 
to a disturbance of a whole arcade of loops. While three distinct 
loops can be identified in Fig. 0(a) near y= 23", 26", and 29", in 
the following we will concentrate on the latter one. This shows 
the clearest oscillating pattern that lasts from ca. £=180 s to 450 s 
and is marked by a black box in Fig. [2(a) and (b). 

To get a clear signal, we average the Doppler shift in the 
oscillating loop in the black box in Fig. [2(b) in the y-direction. 
The resulting (mean) Doppler velocity as a function of time is 
shown by the black symbols in Fig. [2 (c) which clearly shows 
a damped oscillation. To extract the period P of the oscillation 
and its damping time r we fit this oscillation by an exponentially 
damped sinusoidal function. 


fit ) 


Aq exp 


£ — to 


27T 

T 

sin 

-pit-to) 


+ Ait + A 2 . (3) 


Here £ is the time, to the initial time, and Aq the amplitude of the 
exponentially damped sinusoidal function. In addition, A 1 and 
A 2 account for a linear background. Fitting the Doppler shifts 
in Fig. |3] (c) with a uniform weight we obtain a period and a 
damping time of 


P = 52.5 s 

t = 125s 


(4) 
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We applied the same analysis also for some other slit positions. 
When the slit is located midway from the loop footpoint to the 
apex, the average Doppler perturbation in the loop appears very 























F. Chen & H. Peter: Estimate the magnetic field strength in a realistic coronal model 






35 
30 
25 
20 

0 .'100 200 300 400 500 '600 

time [s] 




Fig. 3. Oscillation of the loop. Panels (a) and (b) show the intensity and 
Doppler shift of the Fe XIII line synthesised from the coronal model 
as a function of time and space along the slit. The location of the slit is 
indicated in Fig.[2]and the range of Doppler shifts is from —15 km s“ 1 to 
+ 15 km s _1 , just as in Fig.[2] The black boxes in these space-time plots 
indicate the position of the loop with the clearest oscillation pattern. 
Panel (c) shows the Doppler shift of the loop as a function of time. Here 
the Doppler shifts are averaged along the slit across the loop, i.e. along y 
within the black box in panel (b). The red line shows the fit by a damped 
sinusoidal function as defined in Eq. 0. See Sect. 13.31 


similar to that at the apex, albeit with a reduced amplitude. Still, 
the damping time and period derived at these other positions are 
fully consistent with the results at the loop apex. If the slit is 
positioned close to the loop footpoints, the oscillatory Doppler 
signal gets mixed with slowly changing Doppler shifts due to 
flows (along the mag netic field) filling and draining the loop 
(cf. IChen et al.l2014l) . Therefore, when analysing the oscillation 
close to the footpoints we first remove this smoothly varying sig¬ 
nal before fitting the oscillation by Eq. (j3j. Again the period and 
damping time are consistent with the values found near the apex 
listed in Eq. (Q}. This underlines that the loop shows a coherent 
oscillation all along (see also Sect. l3.4l >. 

Comparing these values of Eq. (jTJi to observations, they are 
found at the lower en d of the distribution foun d in the compila¬ 
tion of oscillations bv iVerwichte et ali i ( 2013b ). This is not sur¬ 
prising, because the the size of our computational box is limited, 
and thus the loops we can study have lengths of below 50 Mm. 
The observed loops are mostly longer by a factor of two (or 
more). Hence the expected period would be longer by that factor 
(the period is proportional to the loop length, see Sect. 14.11 


An important observational test for our model loop is 
the sca ling between period and damping time as derived by 
iVerwichte et al.l (i2013bl) based on a sample of 52 oscillating 
loops. They found r = ctP 7 , with logioa = 0.44±0.31 and 
7 = 0.94±0.12. The period and damping time of the oscillation 
found in our model listed in Eq. (Q} fits this scaling relation very 
well. This suggests that our model for the loop oscillation and its 
damping capture the right physical processes, even though the 
period and damping time are at the lower end of what is found 
in observations. 

3.4. Mode of the oscillation 

Apart from the measurement of the oscillation parameters, the 
identification of the corresponding wave mode is equally impor¬ 
tant. In the synthetic EUV images the loop oscillates primarily 
in the transverse direction without any visible (stationary) node 
along the loop (cf. movie with Fig. [7}. This is supported by the 
Doppler patterns that have the same sign all along the loop (see 
Fig. [2}. Together this suggests that the oscillation in the model is 
a fundamental kink mode. 

A quantitative test for the presence of the fundamental kink 
mode is given through the phase difference of the velocity distur¬ 
bance measured at different positions along the loop (D. Yuan, 
private communication). Because in the fundamental kink mode 
the loop should oscillate coherently all along, this phase differ¬ 
ence should vanish. 

To estimate the phase difference between the velocity dis¬ 
turbance at different positions along the loop we cross correlate 
the oscillation at the apex to four other locations, two midway to 
the footpoints and two close to the footpoints at each side. For 
this we use the fits to the variation of the Doppler shifts from 
Eq. (f3} and show the cross correlation as a function of the time 
lag in Fig. [4] All cross correlations peak at about zero time lag, 
i.e. the velocity disturbances of the oscillation at different posi¬ 
tions have no phase difference 0. Moreover, the peak values of 
the cross correlations are close to unity, indicating that the os¬ 
cillations are very similar at different positions, i.e. that the loop 
oscillates as a whole. 

Another indicator for the fundamental mode is given thor¬ 
ough the variation of the normal velocity along the loop. Because 
for the fundamental mode two nodes of the oscillation should 
be found at the footpoints, there should be a smooth drop of 
the normal velocity from the apex to the feet. This is confirmed 
by the plot in Fig. [5] which shows the normal velocity at a time 
near maximum velocity. There we also illustrate the difficulties 
that would be encountered in real solar observations for this test. 
While in our model we have direct access to the normal veloc¬ 
ity, in observations this is not the case. To show this, in Fig.0 
we also plot the Doppler shift as it changes along the loop when 
observing from straight above (i.e. the shift along the dotted line 
in Fig-EJ- Essentially this shows the vertical velocity in the loop. 
While near the apex this basically reflects the line-of-sight com¬ 
ponent of the normal velocity, near the footpoints this signal is 
heavily contaminated by the flows in or out of the loop. Here 
we see a draining of the loop, i.e. downflows at both footpoints. 
The downflows (12kms _1 and 24 km s 1 ) are comparable with 
the transverse oscillation speed (14 km s -1 ) at the loop apex. In 
a real observations this might prevent any solid conclusions on 


2 For kink mode oscillations, there is a phase difference of ty/2 be¬ 
tween the velocity disturbance and the displacement at the same posi¬ 
tion of the loop. 
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time lag [s] 


Fig. 4. Cross correlations of the oscillation at different locations along 
the loop. The five curves show the cross correlation of the loop apex 
with five positions along the loop (apex, two positions midway from 
apex to loop, and two near the footpoint at both sides). The solid line 
shows the self-correlation at the loop apex. The range of the time lag 
shown here is about four times the oscillation period of 52.5 s at the 
loop apex. The red diamonds indicate the peaks of the cross correlation 
functions. See Sect. 13.41 



normalized arc length 


Fig. 5. Line-of-sight Doppler shift and normal velocity of the oscilla¬ 
tion along the loop. The blue curve shows the velocity of the oscillating 
loop in the normal direction at a snapshot of maximum velocity am¬ 
plitude (at t=200 s). The red line displays the Doppler shifts from the 
synthesised spectral profiles, basically the vertical velocity in the loop, 
which also show signs of the field-aligned flows near the footpoints. The 
arc length has been normalized to the loop length, i.e. the footpoints in 
the photosphere are at 0 and 1, the apex is around 0.5. See Sect. 13.41 


the wave mode by checking the Doppler shift variation along the 
loop. 

This result illustrates also an important constraint on simpli¬ 
fied models for loop oscillations. Because one will never prevent 
such in- or outflows neat the loop footpoints in a complex enough 
model (we never see static loops in the 3D models) or on the real 
Sun, these flow structures have to be taken into account. In par¬ 
ticular, because these flow speeds reach a good fraction of the 
sound speed. 

Together with the the fact that the same periods and damping 
times are found all along the loop (Sect. 13.3k the zero time-lag of 
the oscillation at different places along the loop and the clear ex¬ 
istence of nodes at the footpoints underlines that the oscillation 
in our model loop is indeed a fundamental kink mode. 


Table 1. Parameters of the oscillating loop. 


parameter 

symbol 

value 

from 

coronal loop length 

L 

45 Mm 

EUV 

temperature 

T 

1.5 MK 

model 

internal density 

Hi 

5.7xl0 8 cm -3 

model 

external density ratio 

n e /rii 

0.12 

model 

oscillation period 

P 

52.5 s 

Eq.© 

damping time 

r 

125 s 

Eq.© 

kink mode speed 

Ck 

1730 km/s 

Eq.© 

inverted B\ 

^kink 

79 G 

Eq.© 

average \B\ 

(B) 

92 G 

Eq.© 


Notes. The loop length is derived from the images of the synthesised 
corona (EUV). The temperature, the density in the loop, and the ratio 
of the density outside to the density inside the loop are derived from the 
3D MHD simulation (model). The other quantities are derived from the 
equations referred to in the table. 


4. Oscillation frequency and magnetic field 

4.1. Estimate of the magnetic field strength from coronal 
seismology 

In a seminal naner lEdwin & Roberts! (il983t) showed that the pe¬ 
riod of a loop oscillating in the corona would depend on the mag¬ 
netic field, which in turn allows to estimate the field strength if 
one can identify the oscillation and determine its frequency. In 
an observation the phase speed Ck of a kink mode oscillation usu¬ 
ally is estimated by 

2 L 

Ck = -p (5) 

where L is the loop le ngth and P the oscillation period (e.g 
iNakariakov et al.l 11999 ). The loop length is determined from 
EUV imaging observations and the period by fitting the oscilla¬ 
tion in the displacement of the loop position of the Doppler shift 
(e.g. as described in Sect. 13.3I >. TableQ]summarises these param¬ 
eters derived from our synthesised observation of the oscillating 
loop and other loop parameters and results. The value we de¬ 
rive here for Ck of 1730kms _1 is still below the limiting Alfven 
speed for the simulation, see discussion with Eq. ©. When we 
derive the magnetic field strength, we still consider the limiting 
in the analysis bel ow in Eq. ©. 

As derived by lEdwin & Roberts ( 198 3). the phase speed of 
this kink mode in a slender magnetic flux tube with uniform den¬ 
sity depends on the Alfven speed and the density p inside and 
outside the flux tube, 

/ 2 2 \ 1/2 

(6) 

V P> + Pe > 


Here the subscripts i and e refer to internal and external, i.e. to 
inside and outside the flux tube. 

Because the corona is a low-/? plasma, it is usually assumed 
that the magnetic field strength /? k i n k (derived from the kink 
mode oscillation) inside and outside the loop is the same. In our 
MHD simulation we can confirm that this is indeed the case. 
Using the effective Alfven speed in our model from Eq. ([T]) to 
account for the Alfven speed limiting, Eq. ([6]) reads 


-Bkink 

( 2 f 

/Ai + /Ae 

VPOPi 

V 1 + Pe/Pi J 

2 
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Fig. 6. Plasma parameters and magnetic field strength along the loop. 
Panel (a) shows the temperature (black) and density (blue) for the condi¬ 
tion inside the EUV loop (solid lines labeled internal). The correspond¬ 
ing values just outside the loop are plotted as dashed lines (labeled ex¬ 
ternal). Panel (b) shows the magnetic field strength along the loop in the 
3D MHD model. The red solid line shows the coronal magnetic field 
strength derived from the kink oscillation, Brink, and the red dashed line 
the average magnetic field strength, (B), defined by Eq. {8}. The po¬ 
sition along the loop is normalized by the total loop length (53 Mm). 
The positions 0.0 and 1.0 are at footpoints in the photosphere. The grey 
areas indicate the lower atmosphere, where T< 1 MK. The loop length 
in the coronal part, i.e. between the gray regions, is about 45 Mm. See 
Sect. l3.3l for panel (a) and Sects. 14. 1 I and l4.2l for panel (b). 


The last term [• • • ] is not present when observers analyse their 
data. It is an artifact of the Alfven speed limiting we applied in 
our model (see Sect. [2]) and consequently we have to correct for 
this here, too. 

If observations provided the kink mode wave speed through 
Eq. ([3 and if the internal and external densities can be estimated 
from observations, too, then Eq. ([3 provides an estimate for the 
magnetic field strength in the loop based on the analysis of the 
kink mode, f?ki n k- In absence of a reliable tool for density diag¬ 
nostics when using imaging observations, one often simply as¬ 
sumes a typical value for the coronal density, e.g. 10 9 cm -3 , and 
a typical density contrast of p e /pi = 0.1. 

In our study we use the values for the internal and external 
densities as derived directly from the MHD model (see TableQJ. 
For this we choose a group of magnetic field lines within the 
cross section of the oscillating loop as seen in EUV emission 
for the internal properties and another group of magnetic field 
lines in the ambient corona outside for the external properties. In 
Fig. [6] (a) we plot the number density and temperature averaged 
for the respective group of fieldlines along the loop, separately 
for the conditions inside and outside the loop. 

In the coronal part (i.e. where T>10 6 K) of the EUV loop, as 
well of the flux tube just outside the loop, the temperature and 
the density are roughly constant. This is because of the highly 
efficient heat conduction and the large pressure scale height at 
high temperatures. In this coronal part (between the gray hatched 
areas of Fig. [3, the average density inside the EUV loop is about 
Pi=3.8 x 10~ 12 kgm -3 (nj=5.7 x 10 8 cm -3 ). Outside the loop, 
the external average density is p e =4.4 x 10~ 13 kg m~ 3 (n e =6.5 x 


10 7 cm '). The resulting density ratio is about p e /p ; = 0.12. The 
density we find in the loop is comparable to typically assumed 
values, and also the ratio is similar to the typical assumptions 
applied for tube models and for the inversion of observations 
(e.g. Pe/pi = 0.1 bv lNakariakov & Ofmanl2QOll) . 

The loop length in the coronal part is roughly 45 Mm, which 
is consistent with the length one would derive from the syn¬ 
thesised EUV images. Together with the oscillation period of 
P = 52.5 s derived from fitting the oscillation in Sect. 13.31 we ob¬ 
tain a kink-mode phase speed of Ck= 1730km s” 1 from Eq. ([3 
With the densities derived above we can use this value of Ck 
to derive the coronal magnetic field strength based on the kink 
mode oscillation from Eq. ([3 which yields f?ki n k = 79 G. 

This derivation of the average coronal magnetic field strength 
basically copies the procedure applied to observations. Now the 
main question is what this average actually means. 


4.2. Comparison with the actual magnetic field strength 
along the loop 

The magnetic field strength varies along the coronal loop in both 
the real corona and our numerical model, basically reflecting the 
expansion of the magnetic field. Therefore, the value deduced 
from coronal seismology does not necessarily represent the field 
strength at any particular position of the loop, but is an average 
through the coronal part of the loop. Most importantly, the sig¬ 
nificant variation of the magnetic field along the loop is usually 
neglected in coronal seismology, even though the i mportan ce of 
this effect has been shown before (e.g. lOfman et al.ll2012 1 The 
profile of the magnetic field strength along the loop in our model 
(along the center field line of the loop) is plotted in Fig. [3(b). In 
the the following we will investigate the field strength (e.g. its 
average) only in the coronal part, i.e. where T>10 6 K between 
the gray areas in Fig. [3 

The arithmetic mean of the magnetic field strength along 
the loop is 120 G, i.e. significantly higher than the value of 
/jkink = 79 G derived from the oscillation. The reason for this dis¬ 
crepancy is that coronal seismology estimates the magnetic field 
strength from the observed average wave speed as defined as 
Eq. ©. Because the actual wave speed is (usually) not constant 
through the loop, this gives more weight to regions of low Alfven 
speeds, i.e. to regions where the propagating wave packet dwells 
the longest. These are the regions of low magnetic field strength, 
which is why the magnetic field strength derived from oscilla¬ 
tions is much lower than an arithmetic mean. 

To define an average that is closer to the value derived 
from oscillations, one has to account for the variation of the 
wave speed (or in other words, the wave travel time). Similar 
to lAschwanden & Schriive rl (1201 ll) we define 


(B) = L 



ds 


-l 


B(s) 


( 8 ) 


where s is the coordinate along the arc length of the loop. 
lAschwanden & Schriiveil (i2011! ) argue that a flux tube with a 
field strength B(s) changing along the loop would oscillate 
with the same frequency as a flux tube with a constant field 
strength (B) as derived from Eq. ©. In our model loop, we find 
(B)= 92 G, as indicated by the red dashed line in Fig. [3(b). Com¬ 
pared with the arithmetic average along the loop, this is closer to 
the the magnetic field derived from the kink mode, /Ikmk= 79 G, 
but it is still significantly larger by some 15% to 20%. 

That the average (B) based on Eq. ® comes reasonably 
close to the value derived from the oscillation, /7ki„k, underlines 
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that the magnetic field strength estimated by coronal seismology 
represents the average magnetic field weighted with the wave 
speed along the loop. Basically coronal seismology returns the 
value of the magnetic field that one would find in a flux tube with 
constant magnetic field with the same Alfven crossing time. 

However, the field derived by coronal seismology, /4i,,i^ is 
still 15% to 20% smaller than the average ( B ), i.e. the true mag¬ 
netic field in our simulation (or in the real Sun). This implies 
that the actual wave speed (in the model) is underestimated by 
the theoretical value from coronal seismology. A difference of 
up to 50% between the derived magnetic field strength and the 
actual field strength was also found by iDe Moortel & Pascoel 
( 2009 ) who investigated a MHD model (albeit with a less 
realistic magnetic setup and much simpler thermodynamics). 
They concluded that this is because “the combined effect of 
the loop curvature, the density ratio, and aspect ratio of the 
loo p appears to be mor e im portant than previously expected” 
(IDe Moortel & Pascoell2009h . 

Further more the (magnetic) complexity of coronal loops in 
our model and on the real Sun could lead to such a deviation. One 
possible explanation might be the aspect ratio (i.e. width/length) 
of the loop. The particular oscillating loop highlighted in Fig.Q] 
has an aspect ratio of about 0.04 (2 Mm width, 45 Mm long; 
other loops in the model are thinner). In observations the as¬ 
pect ratio is a bit smaller, typically 0.02 (2 Mm width at 100 Mm 
length). However, it has to be noted that here the width refers to 
the width of the EUV loop, while 3D models show that the cor¬ 
responding density structures these EUV loops are em bedded in 
are thicker dPeter & Bingerti[2012: Chen et al.ll2014l) . Thus we 
can e xpect the aspect r atio to be several percent. 

lEdwin & Roberts (1 9831) noted that the actual kink mode 
speed is equal to Ck defined by Eq. © only if the width of 
the flux tube is much smaller than the wavelength of the dis¬ 
turbance. Otherwi se, it decreases with an increasing aspect ra¬ 
tio (see Eq. 15 in lEdwin & Roberts!! 198.3 ). Therefore it might 
be that in our model (and on the real Sun) this assumption 
of lEdwin & Robertsl ( 1983 ) is not fully applicable. Other ef¬ 
fects, e.g. that plas ma-/3 is non-zero (see Eq. 11 for large j3 
of lEdwin & Robertsil 19831) . a smooth density profile across the 
loop, or flows along the loop might play a role, too, in our model 
as well as on the real Sun. A detailed theoretical investigation 
of these effects is beyond the scope of this paper, but should be 
addressed in a future study. 

5. Damping of the oscillation and Reynolds 
number 

So far we concentrated on the frequency of the loop oscillation 
to investigate the magnetic field strength in the loop. We now 
turn to the damping time r of the oscillation as defined through 
Eq. (© and discuss consequences for the dissipation and thus for 
the (magnetic) Reynolds number. 

Damping is commonly found in transverse loop oscillations. 
The damping time is particularly important, because it reflects 
the rate at which the wave energy is either converted into an¬ 
other wave mode or dis sipated. Using a scaling r elation derived 
from numerical models. iNakariakov et al.1 ( 1999 ) found that the 
magnetic Reynolds number deduced from the observed damping 
time is of the order of 10 6 . This is several orders of magnitude 
smaller than the classical value that would be of the order of 10 14 
and consequently would indicate an anomalously high dissipa¬ 
tion. Later. IQfman & Aschwandenl d2002l) quantified this further 
through an analysis of the scaling relation between damping time 


and loop parameters. They estimated an anomalously large vis¬ 
cosity of 10 9 - 2±J - 5 m 2 s -1 , which also gives a Reynolds number 
that is much smaller than the classical value. These results imply 
that the wave energy could be efficiently converted into heat be¬ 
cause of the anomalously high diffusivity. Still, one has to keep 
in mind, that the damping is not necessarily due to dissipation 
of the wave energy (even though ultimately the wave will be 
dissipated). iGoossens et al.l d2002l) argued that the observed fast 
damping of transverse loop oscillations can be explained by the 
damping of quasi-mode kink oscillations. Thus, they suggested 
that there might be no need for the anomalously large viscosity. 

In our model the magnetic resistivity ry is 5 x 10 9 m 2 s” 1 (and 
the kinematic viscosity v is of the same order). Essentially this 
val ue is determined by the grid spaci ng in the numerical mod¬ 
els dBingert & Peteil201 It lPeterll2015l) and ensures that currents 
that build up due to the driving of the magnetic field by pho- 
tospheric motions will be dissipated at scales of the grid spac¬ 
ing. These values of resistivity and viscosity are in the middle of 
the range of values found bv lOfman & Aschwanden ( 120021) . i.e. 
10 9 '2±3.5 m 2 s -i p rom this we conclude that the dissipation co¬ 
efficient in our model (10 9 7 m 2 s -1 ) is consistent with the high 
values found in observations. 

If we take the loop half width of our model loop («1 Mm) 
as a length scale, the dissipation time for the wave in our model 
would be Tdi ss ~ (1 Mm) 2 /5x 10 9 m 2 s -1 = 200 s, which is close 
to the damping time of 125 s for the oscillation in our model. 
This suggests that in our model the damping is indeed due to 
(resistive and viscous) dissipation. A further more detailed anal¬ 
ysis of the energy budget will be needed to investigate how much 
of the energy released by the trigger is first converted into the os¬ 
cillation and then dissipated by viscosity and resistivity. In par¬ 
ticular further investigations will have to study the role of the 
spatial resolution and thus of the Reynolds number on the result¬ 
ing oscillations. 

6. Conclusions 

We presented a loop oscillation found in a realistic coronal 
model driven by magnetic flux emergence through the photo¬ 
sphere (i.e. the bottom boundary). The treatment of the energy 
balance in the model allows us to synthesise imaging and spec¬ 
troscopic observations from the model. There we found a clear 
transverse loop oscillation that we identified as the fundamen¬ 
tal fast kink mode. The period (P= 52.5 s) and damping time 
(t=125s) of the oscillation are consistent with observations. 
In particular, our oscillation also follows the o bser ved relation 
between period and damping time deduced by IVerwichte et al.l 
(l2013bl) . At least in our model the damping of the oscillation 
is due to resistive and viscous dissipation. To what extent this 
conclusion can be also drawn for the real Sun will have to be 
revealed by new simulations at higher spatial resolution. 

We applied methods of coronal seismology to the coronal 
emission synthesised from our model corona in the same way as 
it would be done for real observations. Here we concentrated on 
the oscillation in the Doppler shift of a coronal emission line. We 
chose a (vertical) line-of-sight for the analysis of the model data 
that mimics the observation of a coronal loop near disk center. 

Based on the fundamental kink mode oscillation we deduced 
an average magnetic field strength of f?kink=79G in the loop. In 
contrast to solar observations, in the model we know the mag¬ 
netic field in the oscillating loop. This way we can understand 
what the field strength derived from coronal seismology actu¬ 
ally means. In the coronal part of the model loop the magnetic 
field strength varies strongly by about a factor of five and drops 
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to some 50 G near the apex, i.e. significantly below the value 
derived from coronal seismolog y, /j k ink - Th e ave rage of the field 
along the loops as suggested bv lAschwanden & Schriiverl (1201 ll) 
gives a relatively good match, (B)= 92 G. This average of the 
field corresponds to the constant field in a flux tube with a con¬ 
stant cross section that has the same wave travel time. Still, the 
difference between /j knlk and (B) is considerable. Because many 
of the assumptions of the simple derivation of the wave period 
are violated in the model (and most certainly on the real Sun), 
e.g. the loop width, substructure, or dynamics, this is not too 
surprising. 

We conclude that the magnetic field strength deduced by 
coronal seismology can be a good representative of that in the 
upper part of coronal loop. Realistic models like the one pre¬ 
sented here might further guide the interpretation of coronal os¬ 
cillation results, in particular to relax some of the assumptions 
to get a better match of coronal properties derived from coronal 
seismology. 
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